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We study the effect of monomer immobihzation (quenching) on the orientation of the lamellae 
in symmetric diblock copolymer thin films with neutrally wetting surfaces. A small fraction of 
the monomers immediately next to the solid substrate is presumed to be quenched. In both the 
weak segregation limit and the strong segregation limit, quenching favors the lamellae orienting 
perpendicular to the film. Quenching inhibits the order-disorder transition twice as much for the 
parallel orientation as for the perpendicular. 
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0> ! I. INTRODUCTION 

^ ■ Bulk diblock copolymer melts have been studied extensively, both experimentally and theoretically!]'! Recently, 

Qmany studies have examined thin films of diblock copolymers. In this paper, we focus on symmetric diblock copoly- 
mers. In the bulk, symmetric diblock copolymers form lamellae upon microphase separation. In thin films, a varijetK 
lO of behavior has been observed. In some instances, the copolymers form lamellae parallel to the film interfaces. lTEI 
(See figure 1(a).) Each of the two interfaces usually preferentially attracts one of the copolymer species. Having the 
" ' ■ lamellae parallel to the integfaees allows the maximum amount of favorable contact. But in many experiments, other 
, morphologies_a,ce. observed .Q^liil (See, for example, figures 1(b) and 1(c).) Some of the results can be explained by 
confinement .E3ilj Frustration resulting from the film thickness being incommensurate with the natural lamellar period 
I can cause the morphology of figure 1(a) to become energetically unfavorable. However, confinement does not appear 
fS| • to be sufficient to explain all the experimental observations. In this paper, we study the effect of monomer mobility 
T-H on the morphology. There is reason to believe that, at least _lh some situations, the mobility of some monomers 
, immediately adjacent to a solid substrate is greatly reducedEOilj To make the calculation tractable, we assume that 
' such monomers are completely immobilized (that is, quenched). 

Our system is a thin polymer film of thickness h. Each polymer has N/2 monomers of some species A and N/2 
monomers of species B. The z = surface of the film is adjacent to a solid substrate, while the other surface z = h 
I is free (in contact with air, for example). Since we are interested in the effect of quenching, not confinement, we 
. assume that the film can adjust its thickness to eliminate any incompatibility between the film thickness and the 
(~| natural lamellar period. We assume that the film thickness h is moderate — large compared to the chain root-mean- 
O , square end-to-end distance Re but not so large that the bulk properties overpower the interfacial effects. A small 
O • fraction of monomers immediately adjacent to the solid substrate at z = is assumed to be quenched. To simplify 
^ . the calculation, we assume that the fraction of quenched monomers is sufficiently small that at most one monomer 
per chain is quenched. We also assume that it is equally likely for A and B monomers to be quenched. Finally, we 
assume that both the solid substrate and the free surface are neutral — that is, there is no preferential segregation 
^ of A or B monomers to the substrate or the free surface. This last assumption differs fmiHJ.ie conditions in most 
■ ■ ' experiments, though through careful design, it is possible to satisfy this last assumptionU'EHltiJ In future studies, we 
would like to eliminate this final assumption. We consider both the weak segregation limit and the strong segregation 
limit. Previous theoriesE^I have treated the effect of surfaces on microphase separation, but they have not considered 
the effect of quenching. We show that monomer immobilization alone is sufficient to change the phase transition. 
It should be noted that the ability of quepchiiig to alter the properties of materials in nonobvious ways has been 
demonstrated in a large variety of systems. Ej~E£l 
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II. WEAK SEGREGATION LIMIT 

Our fiormalism for the weak segregation limit is that of Marko and WittenE^'B and is similar in spirit to that of 
Leibler.LJ First, consider a reference (unquenched) system in which the A and B monomers are chemically identical 
— that is, each polymer chain consists of N identical monomers, the first half labelled A and the other half labelled 
B. The partition function for this system is 

Z,ef = / dRe-^[^l (I) 
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where J dR. is an integral over monomer positions of all chains and ^[R] is the free energy arising from chain 
connectivity. The equilibrium average of an arbitrary quantity X[R] is given by 

/dRe-^Wx[R] 
- JdRe-^W 

Next, apply a perturbation free energy that arises from the immiscibility of A and B monomers: 

5' - - / drM(r)p(r) (3) 



where r = {x,y,z) denotes the three spatial coordinates, p{r) — pA(r) — Pb{y) is the order parameter for phase 
separation, and /9A(r) and pb{^) are the local volume fractions of A and B monomers respectively. (Note that 
/9A(r) + pB(r) = 1.) The partition function for the perturbed system is 



Z = y dRe-^[^le^' (4) 

and the equilibrium average is given by 

,^,_ /dRe-^We^'X[R] 
- /dRe-^We^' 

Since the weak segregation limit corresponds to small perturbations (small /i), we may express the order parameter 
as an expansion in the external potential /i: 

(p(ri)) = ||^-G«(ri)+ J dr2G(2)(ri,r2)/i(r2) + 0(/i') (6) 
where G'(^^(ri) and G*^^^(ri,r2) are defined as 



(9^(ri) 
and 

InZ 



/x=0 



G(2)(ri,r2) = 



(7) 



(8) 



a/i(ri)a/i(r2 
Using equation ^, we find that 

G«(ri) = (p(ri)),.ef-0 (9) 

and 

G(2)(ri,r2) = (p(ri)p(r2))..cf - (p(ri)),.cf (p(r2))rcf - (p(ri)p(r2))rof (10) 
Thus equation ^ can be rewritten (neglecting fi^ and higher order terms) as 

(p(ri)} = J dr2G(^)(ri,r2)/x(r2) (11) 

In our copolymer melt, we use the Flory-Huggins form for the free energy perturbation!^ 

5PH - A/drpA(r)pB(r) ^ - ^ / drp^ (r) + constant (12) 

The demixing parametei0 A has units of inverse volume and is related to the Flory x parameteiS by AV^ — xiV, 
where V is the chain volume and N is the number of monomers per chain. Using a mean-field approximation, 
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and equation [Til becomes an integral eigenvalue equation: 

(P(ri)) = ^ J rfr2G(2)(ri,r2)(p(r2)) (14) 

For small values of A, the only solution to equation |lj is (p(r)) — everywhere — that is, (pA(r)) = (PB(r)) 
everywhere, and the copolymer melt is structureless. The smallest value of A which yields a nontrivial (p(r)) marks 
the onset of microphase separation. 

For a bulk copolymer melt, the results are well known. G^^^(ri,r2) = 2(pA(ri)/CA(r2)}rcf — 2(pA(ri)pB(r2)}rof is 
readily calculated from random walk statistics. There can be no contribution to the (• ■ •)'s unless ri and r2 are on the 
same polymer. Thus the averages can be calculated by summing over all possible pairs of monomers along a single 
chain. It is convenient to label a monomer by the volume v displaced by the section of chain between that monomer 
and the A end of the chain. We label the monomer at ri by vi and the monomer at r2 by V2. Thus 

9 l'V/2 i-V/2 / \ 3/2 

V Ja Jo \2tt\vi-V2\J 

fV / „ \3/2 

dvj dvJ— : e-^^l"-^-"-^! (15) 

V Jo Jv/2 \2n\vi-V2\J 

where a is the "packing length"!! and V is the volume displaced by an entire chain. Since in this case G(^^(ri,r2) 
depends on ri and r2 through ri — r2 only, the eigenvalue equation |lj can be readily solved using Fourier transforms 
and the convolution theorem: 

(p(k)) = ^G(2)(k)(p(k)) (16) 

The onset of microphase separation occurs at AV = 10.495 and k* = 4.77/i?e, where Re is the root-mean-square 
end-to-end distance; (p(r)) can be any sinusoidal function with a wavevector of magnitude k* . A convenient basis set 
for {p{r)) consists of functions of the form 

[kx >0]] ( cos{kyy) [ky > 0] ] ( cos{k^z) [k^ > 0] 
M-r) = 1 [kx^O] } x{ 1 [ky^O] \ X I 1 [k,^0] } (17) 

[kx > 0] J [ sin(feyy) [ky > 0] J [ sin(A:,z) [k, > 0] 

where the factors in the three braces are multiplied together in all possible combinations such that k'^ + ky + k^ is 
satisfied. (There are 26 total combinations; 1x1x1 does not satisfy k'^ + ky + k'^ = fc*^.) Thus (p(r)) — J2i <^i''l'i{^)i 
where the Ui are arbitrary constants. To obtain further information about the microphase separation, we must consider 
higher order terms not included in our formalism. LeiblernJ has calculated these higher order terms for all the possible 
morphologies consistent with k"^ + ky + k^ = k*^ and has found that the lamellar morphology is the most stable. 

In the calculation described below, we calculate the effect of quenching before taking into account Leibler's higher 
order terms. This is reasonable if the perturbation to the microphase separation due to quenching is larger than the 
perturbation resulting from the higher order terms. We can also apply the perturbations in the opposite order. That 
is, we can first consider the higher order terms, which tell us that the most stable morphology is lamellar. We can 
then compare the effect of quenching on the different lamellar orientations. Such a calculation would be very similar 
to to the one described below, and the final result does not change. 

In a thin film, the boundaries at z = and z = h modify the bulk-caLculation slightly. We use Silberberg's 
chain-swapping procedure, which leads to reflective boundary conditions Ho Thus each of the two Gaussian terms 




in equation uB should be replaced by exp 



I |2 
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where ri is the reflection 



21-01 — f 2 I 

of ri in the z = plane. Also, as noted above, the film thickness h is free to adjust. With these two assumptions, we 
find that, as in the bulk melt, AV = 10.495 and k* = A.77/Re. However, the reflective boundary condition at z = 
restricts the basis set functions V'i(r) of equation ^ to functions that are even with respect to z = 0; that is, in the 
thin film, (p(r)) is a linear combination of 

[ sin(fc,x) [kx > 0] J [ sin(fc^y) [fc^ > 0] J ^ ' ~ J 

where k'^ + ky + k^ — k*^. The reflective boundary condition at z = ft. restricts ilJi{r) to functions that are even with 
respect to z — h] thus kzh — ni:, where n is a nonnegative integer. 
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The discussion above applies to unquenched systems. Applying the same formalism to quenched systems results in 
modifications to equations]!^ and ^ 



(r 1,1-2) = (p(ri)p(r2))rcf,q - (p(ri))rcf,q(p(r2))rcf,q (19) 

(^M; = ^ 1 dr^Wij^W^^ (20) 

where (• • •)ref,q and (• • •)q denote equilibrium averages subject to the constraint q and the overbar denotes averaging 
over the constraints. In our system, the constraint is that a small fraction of the monomers immediately adjacent to 
the solid substrate at z = can not move. 

The first term of Gq^'' (ri , r2) in equation [I^ is equal to the first term of G*^^-' (ri , r2) in equation ^ since for any 
arbitrary quantity X[R], {X),-c{ — (X)rcf,q. However, while the second, residual term of G(^)(ri,r2) is zero, the 

residual term of Gq (ri,r2) is nonzero. Using random walk statistics, we find that for a thin film with quenched 
monomers this residual term is 



(p(ri))rcf,q(p(r2))rof,q 2 (pA (ri ))ref ,q(/OA (l"2 ))ref ,q - 2(pA(ri))ref ,q(pB (r2))rof ,q (21) 



where 



{pA (ri ))ref ,q (pa {^2 ) )ref ,q 



V rV/2 / \ 3/2 



V/2 / \ 3/2 

dvJ - ] ^-^^[i..-.r-+iv2-v'f+4] 



4cr 



2t: \ v — V2 \ ^ 

V pV/2 nVl2 



[(3:1-2:2) +(y\-yi) ] 



e 2|„-„ii^i . \ e 2i„-„2i"2 (22) 



y Jo Jo Jo 27r(|w - wi I + |i; - 1^2 1) 

1/2 , / „ \ 1/2 



2Tr\v~vi\J \2t:\v — V2\ 



and (pA(ri))rof,q(pB(i"2)}rof,q is the same except that the integral of V2 goes from V/2 to V instead of to V/2. The 
quantity a is the number of quenched monomers per unit area at the z = boundary. As noted above, we assume 
that the amount of quenching is sufficiently small that at most one monomer per chain is quenched. The integral 
a J dx' J dy' accounts for the fact that it is equally likely for any boundary location (a;',y',0) to have a quenched 

monomer. The integral y dv accounts for the fact that it is equally likely for any monomer of each chain to be 
quenched. 

Since the quenching density a is small, (p(ri))rpfiq(p(r2)}rcf.q ^ {pi^i)£(l2))rci,q, and we can solve the eigenvalue 
equation ^ using degenerate perturbation theoryO We rewrite equation ^ as 



dr^ 



(p(ri)p(r2)}rof,q - (p(ri))rof,q(p(r2))rof,q (p(r2))q = A(/3(ri)}q (23) 

.(0) 

{0) 



where A = 2/A. In the absence of quenching, the solution to equation 23 is A'^*'^ = 2/A* and {p{r))q ~ J2i fliV'i(r) 
where the of equation |l^ form a basis set for the unperturbed solutions and the ai are arbitrary constants. In 

the presence of a small amount of quenching, we can perform an expansion in a about A^''^ and (p(r))q 

A = A(°)+fTA«+a2A(2) + ... (24) 

-(0) , -r-TT^C^) , 2T^T^(2) 



(p(r))q = (p(r))q^ + a{p{v))^' + a'{p{v))^' ' + ■■■ (25) 

Applying degenerate perturbation theory and requiring that / drtpi{r)ipj {r) — for i ^ j, we can calculate trA*^^^ by 
solving 
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- ^fli J dri y dr2(p(ri))ref,q(p(r2))ref,qV'j(ri)'0i(r2) - aX^^^aj J drip'j{r) = 



(26) 



From equations 18, |T], and we see that J dridr2(p(ri)}icf,q(p(r2)}icf,q'0j(ri)V'i(''2) = for i 7^ jS Thus our 
original eigenlunctions ipi are already adequate to treat the perturbed system, and equation En reduces to 



/ dri / dr2(j0(rl))l.cf,q(/0(r2)}ref,qV'^(rl)V'i(r2) 

Jdrt/jfir) 



(27) 



The integrals / dr, J dr^, and / dr2 range over the volume of the thin film. We assume that the extent of the film 
in the x and y directions is large enough that the x and y boundaries have a negligible effect. In the z direction, the 
integrals run from to h. As noted above, h is assumed to be large compared to the polymer rms end-to-end distance 

Re- Thus it is a good approximation to replace /q dzi and J^^ dz2 in the numerator by dzi and dz2. 
For eigenlunctions with nonzero fc^, 

cos(fc2;a;) [kx > 0] 

M^) = <( 1 [fcx = 0] 

sin(fc2:a;) [kx > 0] 
where + ky + k^ = k*^, wc find that 



cos{kyy) [ky > 0] 

1 [ky - 0] 

sm{kyy) [ky > 0] 



X cos(fc2z) [k^ > 0] 



(28) 



(29) 



and 



A(") +crA(i) 



A(o) 1 



A(o) 



aA(i) = 1 + 



2caT/^A(") 
/J 



(30) 



/hare A*^"-* corresponds to the microphase separation transition in the absence of quenching and c is a constant defined 



1 
1 

0.0204 



V 



V/2 



dv 







dvi 

v/2 



V/2 



dv2e 



dvi 



On the other hand, for eigenfunctions with k^ = 0, 

cos(A:j.a;) [k^ > 0] 
= <( 1 [kx = 0] 

sin(fc^a;) [k^ > 0] 



dv 
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v/2 



cos{kyy) [ky > 0] 

1 [ky = 0] 

sin(fcyy) [ky > 0] 



X 1 [fc, = 0] 



(31) 



(32) 



where k'i 



and 



fc2 



fc? = fc*^, we find that 



A « A(") 1 + 



2cctV^ 



(33) 



(34) 



We are interested in the eigenfunctions with the smallest value of A. Thus, in the quenched system, the microphase 
separation transition occurs at A given by equation ^ and (p(r)) = aiipiir), where the V'i(r) g-cjc given by equa- 
tion [3^ and the are constants. Taking into account the higher order terms calculated by LeibleiO further restricts 
(p(r)) to lamellar morphologies. The only lamellar (p(r)) that can be constructed using the V'i(r) of equation ^ as 
basis functions are (p(r)) = cos(A;a;a; + kyy + 9) where kl + ky = k*^ and is a phase constant — that is, the lamellae 
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are perpendicular to the substrate, as in figure 1(b). (Notice that having lamellae parallel to the substrate — that is, 
(/9(r)} = cos(fc*2:) — would require ■(/'i(r) given by equation 28, corresponding to A given by equation Since this 
is not the minimum value of A, the tpi{r) given by equation 28 cannot be used.) 

Thus quenching inhibits phase separation since quenching increases the value of the demixing parameter A at the 
order-disorder transition. Furthermore, quenching inhibits the formation of lamellae parallel to the substrate twice as 
much as lamellae perpendicular to the substrate. Consequently, quenching favors the lamellae orienting perpendicular 
to the substrate. From equation ^ we see that the amount of inhibition increases with increasing quenching density 
a and decreases with increasing film thickness h. Since quenching is an interfacial effect, increasing the volume of 
material relative to the interfacial area decreases the effect of quenching. We can also estimate the order of magnitude 
of the increase in A due to quenching. A large change in A would require large a and small h. The largest monomer 
quenching density consistent with the assumptions in our calculation is cr « cl/V, corresponding to about one quenched 
monomer per chain touching the substrate. Our formalism can handle higher quenching densities, but equation 
would need to be modified to allow for more than one quenched monomer per chain. The smallest film thickness 
allowed by our calculation is ft. « i?e ~ \/V/a. Having h smaller than Re would cause our calculation to break down 
in several places. Thus the largest possible perturbation consistent with our assumptions is 



A(o) h yi/2 

For a typical polymer, the chain volume V « lO^A and the packing length a w lOA. We also know that c w 0.02 and 
A(°)y « 10, so AA/A^"^ « 0.02. A 2% change in the demixing parameter A could be difficult to detect experimentally 
but is not so small that one can ignore it completely in interpreting experimental results. 



III. STRONG SEGREGATION LIMIT 



In the strong segregation limit, we assume that symmetric diblock copolymers form lamellae consisting of regions of 
A monomers separated by sharp interfaces from regions of B monomers. In a thin film, the lamellae can orient either 
parallel or perpendicular to the substrate. In the absence of quenching, the free energies of these two orientations are 
equal, but this degeneracy is broken when there is quenching. Consider the effect of quenched monomers when the 
lamellae are parallel to the substrate, as in figure 1(a). Half of the quenched monomers are stuck in an unfavorable 
region. In order for a polymer chain containing an unfavorably quenched monomer to reach a favorable region, the 
disfavored chain segment must stretch a distance L/4, where L is the lamellar period, and this costs substantial free 
energy. Now consider the effect of quenched monomers when the lamellae are perpendicular to the substrate, as in 
figure 1(b). Again, half of the quenched monomers are stuck in an unfavorable region. In order for a polymer chain 
containing an unfavorably quenched monomer to reach a favorable region, the disfavored chain segment must stretch 
a distance of at most L/A. Most chains containing an unfavorably quenched monomer need only stretch a shorter 
distance and thus incur a smaller free energy penalty. Therefore, as in the weak segregation limit, quenching causes 
the lamellae to orient perpendicular to the substrate. 



IV. DISCUSSION 



The calculation reported above describes the effect of immobilized monomers on block copolymer phase separa- 
tion. We have shown that this immobilization tends to favor perpendicular lamellae. We have also shown how the 
immobilization constraint can be incorporated into the standard field theory used to treat polymer phase separation. 
On the other hand, our theory does not yet give an adequate account of why perpendicular lamellae are observed 
experimentally. Below we discuss the various limitations in our theory, their experimental relevance, and the prospects 
for improving the theory. 

The most unrealistic assumption is that there is no preferential interaction of the A or B monomers with the 
surface^ jSiich neutral surfaces are uncommon in practice. The reports of perpendicular lamellae that motivated this 
studyLltl'Ej were done on surfaces that were not at all neutral. Preferential interaction with the surface naturally 
favors a morphology in which the surface is covered with the preferred species. It thus opposes the tendency of the 
quenched chains to make perpendicular lamellae. The preferential surface interaction amounts to a perturbation in 
the free energy that is linear in the order parameter p. The perturbation due to quenching, on the other hand, is an 
effect quadratic in p: it does not break the A-B symmetry. Thus, for any significant phase separation amplitude, we 
expect the linear preferential adsorption term to dominate, unless its coefficient is zero (i.e., neutral surfaces). 
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In any case, even if the surfaces are nearly neutral, there are na lack of effects favoring perpendicular orientation, 
even without quenching. The effect of incommensurate thicknesal3~ll3 was already mentioned in theJntroduction. In 
addition, nematic interactions between the monomers and the surface favor perpendicular lamellae.ES Given all these 
restrictions, we do not expect quenching to play an important role in the weak segregation regime. In the strong 
segregation regime, the quenched chains carry a large free energy penalty and could play a much greater role. Our 
analysis of this regime is in progress. 

Another restriction of our theory is the assumption of a small number of immobilized monomers. Under this 
assumption the effect on the phase separation threshold A is of course small as well. Still, the calculation can give a 
glimpse of the effects of larger quenching density. If the quenched monomer density approaches one per chain touching 
the surface, the volume fraction of quenched chains near the surface becomes of order unity. Beyond this point our 
calculation becomes unreliable. As the number of immobilized monomers increases beyond this point, the surface 
chains become quenched at multiple sites, thus forming a series of grafted loops. One anticipates that these loops 
would inhibit microphase separation even more than the singly-attached chains do, since the fluctuations that allow 
phase separation are more inhibited. 

Our results are amenable to quantitative tests. Such tests would be of interest, as they would reveal how immobiliza- 
tion affects phase separation morphology. Naturally, one possible test is to induce block copolymer phase separation 
on a neutral surface containing a few grafted chains. Though we have only treated chains immobilized at an arbitrary 
monomer, our method is easily adapted to treat various cases — for example, end-grafted chains or chains grafted at 
the A-B junction point. One need only replace the average y dv in equation by the appropriate superposition 
of V values. An experimental test could reveal the optimal way of influencing the morphology. It could also extend the 
understanding of the quenching effect beyond the narrow limits treated above. Such an experimental test, however, 
would face obstacles. Any residual non-neutrality of the surfaces would inhibit perpendicular orientation. In addition, 
nematic interactions with the surface would favor perpendicular orientation even in the absence of quenching, as noted 
above. A cleaner investigation could be performed using computer simulations. Here the neutral surface and the A-B 
symmetry of the chains could be realized exactly. Such simulations could help elucidate the interplay between the 
various factors, including preferential surface interactions, nematic interactions, film confinement, and quenching, that 
influence the phase separation morphology. 



V. CONCLUSIONS 



We have shown that monomer mobility can affect the morphology of diblock copolymer thin films. In our system, 
immobilizing a fraction of the monomers next to the substrate favors the lamellae orienting perpendicular to the 
substrate. In the future we would like to study the effect of quenching in other systems since we believe that 
quenching occurs to some extent in a wide variety of systems. 
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FIGURE CAPTIONS 



Figure 1. Some possible morphologies for thin films of symmetric diblock copolymers, as suggested by experiments. 
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Figure 1 (a) 
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Figure 1 (b) 




Figure 1 (c) 
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